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ABSTRACT 



> , 

In a recent paper Loveday et al. (1992) have presented new results on the luminosity 

^sD . 

function for a sample of galaxies with bj < 17.15. After having morphologically classified 

m 1 

each galaxy (early-type, late-type, merged or uncertain), they have estimated the parame- 
OS ■ ters of a Schechter luminosity function for early- and late-type galaxies. However, in their 



sample there is a bias against identifying early-type galaxies at large distances and/or faint 



magnitudes: in fact, many of the early-type galaxies at faint magnitudes have probably been 
classified as "uncertain". As discussed in Loveday et al., the existence of such a bias is 



indicated by the fact that for these galaxies < V/V max >= 0.32. 

In this paper we show, both theoretically and through the use of simulated samples, 
that this incompleteness strongly biases the derived parameters of the luminosity function 
for early-type galaxies. If no correction for such incompleteness is applied to the data (as 
done by Loveday et al.), one obtains a flatter slope a and a brighter M* with respect to the 
real parameters. 
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1. Introduction 



An accurate knowledge of the optical luminosity function of galaxies is required for many 
applications in cosmology. For instance, it is essential in interpreting galaxy number counts 
and in analyzing the spatial distribution of galaxies from redshift surveys; in addition, the 
shape of this function is of theoretical interest as it may provide constraints on models of 
galaxy formation. 

An interesting question about the luminosity function of galaxies concerns its universal- 
ity: indeed, Binggeli, Sandage & Tammann (1988) have shown that the luminosity function 
depends on the morphological type, particularly at the faint end. On the other hand, it has 
been demonstrated that the mix of morphological types is closely related to the local matter 
density (Dressier 1980). The accurate knowledge of the luminosity function for each mor- 
phological type is of great interest also for the models of number-magnitude counts. These 
models strongly depend on the morphological mix and therefore need the knowledge not only 
of the fraction of the various galaxy types but also of their K-corrections and luminosity 
functions. 

Loveday et al. (1992) have recently presented new results on the luminosity function 
for a sample of galaxies with bj < 17.15. After having morphologically classified each 
galaxy (early-type, late-type, merged or uncertain), they have estimated the parameters 
of a Schechter luminosity function for early- and late-type galaxies, using the STY para- 
metric maximum likelihood method (Sandage, Tammann & Yahil 1979). While for late-type 
galaxies their parameters are in reasonable agreement with those derived from other samples 
(see f.i. Efstathiou, Ellis & Peterson 1988), the parameters for early-type galaxies are not 
consistent with previous determinations. 

As mentioned by Loveday et al., in their sample there is a bias against identifying 
early-type galaxies at large distances and/or faint magnitudes: in fact, many of the early- 
type galaxies at faint magnitudes have probably been classified as "uncertain" , and therefore 
have not been used in computing the luminosity function. The existence of such a bias is 
demonstrated by the fact that for these galaxies < V/V max >= 0.32. The same bias appears 
not to be present in the classification of the late-type galaxies, for which < V/V max >= 0.47 
(see Table 1 in Loveday et al. 1992). In this paper we show, both theoretically and through 
the use of simulated samples, that this incompleteness strongly biases the derived parameters 
of the luminosity function for early-type galaxies. 
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In Sect. 2 we demonstrate that the classification incompleteness biases the results of the 
STY method and in Sect. 3 we estimate the amount of this bias through simulations. 



2. The luminosity function of galaxies 

The luminosity function of galaxies is well represented by a Schechter (1976) form 



where a and L* are parameters referring to the shape of the function and 0* contains the 
information about the normalization; these parameters have to be determined from the data. 

Many different methods have been used in the past years to compute the parameters 
of the galaxy luminosity function. Recently, however, the STY method (Sandage et al. 
1979) has been the most widely used, and it has been shown that this estimator is unbiased 
with respect to density inhomogeneities (see f.i. Efstathiou et al. 1988). The basic idea 
of this method is to compute the estimator of the quantity ^, where $ is the integrated 
luminosity function. Under the assumption that 4>(L) is not a function of position [i.e. 
4>(L,x,y, z) dL dV = p(x,y, z)dV ip(L)dL], the probability of seeing a galaxy of luminosity 
Li at redshift is 

ft = (2) 



/•OO 

/ ip(L)dL 



in(Zi) 

where L min (zi) is the minimum luminosity observable at redshift Zi in a magnitude-limited 
sample. 

The best parameters a and L* of the luminosity function are then determined by max- 
imizing the likelihood function C{a,L*), which is the product over all the galaxies of the 
individual probabilities pi. This corresponds to minimize the function 
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where V is the incomplete Euler gamma function and iV is the total number of galaxies in 
the sample. 

This formula is correct only for a complete, unbiased sample in which all galaxies with 
Tfi < raum are members of the sample or all galaxies with m < mu m have the same probability 



of being members of the sample (as, for example, in a redshift survey with — sampling). In 

n 

other cases in which each galaxy of the sample has a different weight Wi, which may be a 
function of an intrinsic property of the galaxy (f.i. the distance, the absolute or apparent 
magnitude, the diameter, etc.), eq. (2) is not valid anymore. If we define the weight Wi as 
the inverse of the probability that the i th galaxy has of being included in the sample, eqs. 
(2) and (3) have to be modified in: 

/ \ 
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and 
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(5) 

Loveday et al. (1992) have used eq. (2) to compute the luminosity function for the galaxies 
in their sample, both when they considered all galaxies and galaxies divided in sub-groups, 
as a function of the morphological type. Since, as mentioned in the Introduction, their 
morphological classification of early-type galaxies is biased at faint apparent magnitudes, a 
weight Wi(rrii) should be associated to each galaxy and the use of eq. (2) for determining the 
luminosity function of early-type galaxies is not correct anymore. In the following section 
we will quantify, through simulated samples, the differences between the results obtained 
through the use of eq. (3) and (5). 

3. Results 

In order to estimate in a quantitative way the error introduced applying eq. (3) instead 
of eq. (5) to an incomplete sample, we have used two types of random simulations. 

1) We have randomly distributed 5 millions of points following a Schechter luminosity 
function with parameters a = —1.10 and M* = —19.50 (which are typical values for the 
galaxy luminosity function), obtaining a sample of ~ 11000 galaxies with mu m = 17.15. Then 
we have introduced in this sample an incompleteness function, depending on the apparent 
magnitude of the galaxies in the form of 



f(m) = a(m - m ) 



(8) 
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and for each m we have randomly eliminated from the sample a fraction /(to) of galaxies. 
We assumed that /(to) = for m < m , i.e. for galaxies brighter than m the sample is 
complete. We have chosen m = 16, assuming that for galaxies brighter than 16 th magni- 
tude the morphological classification is relatively easy, and we have used various values for 
a, corresponding to different values of < V/V max > and, therefore, to different levels of in- 
completeness. Then we have computed the parameters of the luminosity function for these 
samples using both eq. (3) and eq. (5). 

Table 1 lists the derived parameters for three representative cases: column (1) gives the 
adopted incompleteness function, column (2) the < V/V max > of the sample and column (3) 
the number of galaxies; the parameters a and M* derived from eq. (3) and (5) are listed in 
columns (4) and (5), (6) and (7), respectively. From this table it is clear that, as a increases, 
the use of eq. (3) produces flatter slopes and brighter M* with respect to the real parameters. 
In all these cases the derived parameters are not compatible with the real ones, as shown 
by the confidence ellipses in Fig. 1. But if we use eq. (5), which takes into account the 
incompleteness function, we obtain "corrected" parameters in very good agreement with the 
real ones. 

2) In the previous case we have used a simulation with a large number of points in order 
to minimize the effects of statistical fluctuations in the test of eq. (3) and (5). Now we try to 
reproduce as much as possible the characteristics of the early-type galaxy sample of Loveday 
et al. (< V/Vmax >= 0.32). We have generated 100 random catalogues, distributed following 
a Schechter function with the parameters derived by Efstathiou et al. (1988) for early-type 
galaxies (a = —0.48 and M* = —19.37). Note that these values are not consistent with the 
values a = +0.2 and M* = —19.71 found by Loveday et al. (see below in Fig. 2). Then we 
have applied the incompleteness function /(to) = 0.8(m — 16), which gives < V/V max >= 0.33 
(this value has been computed as the mean of the values of the < V/V max > derived for each 
catalogue). The number of object in each incomplete catalogue is of the same order as the 
number of galaxies in the Loveday et al. sample (311 objects). Finally, we have computed a 
and M* applying both eq. (3) and (5). The results are given in Table 2, whose columns have 
the same meaning as in Table 1, except that in this case all the parameters are the mean of 
the values derived for each catalogue. In Fig. 2 we show in the M* — a plane the parameter 
pairs found for each incomplete catalogue: open and solid circles refer to parameters derived 
with eq. (3) and (5) respectively, while the star indicates the location of the input parameters 
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and the cross represents the Loveday et al. parameters. The dashed and dotted curves are 
the la confidence ellipses of the parameters as derived by Efstathiou et al. and Loveday et 
al., respectively. The figure shows that there is a clear separation between the two sets of 
points: while the solid circles are very well consistent with the input parameters (star), the 
open circles are displaced toward a flatter slope a and a brighter M*. 

In conclusion, our simulations suggest that, qualitatively, the correction for incomplete- 
ness should move the (a, M*) parameters of Loveday et al. toward those determined by 
Efstathiou et al.. Taken at face value, the shift between uncorrected and corrected values 
as resulting from the simulations is only about half of what is needed to obtain a perfect 
agreement between the results for the two samples. With such a shift, however, the two con- 
fidence ellipses would overlap somewhat at the lcr level, thus becoming reasonably consistent 
with each other. Moreover, it is important to stress that in our simulations we were forced 
to use arbitrary incompleteness function. A detailed knowledge of the functional form of 
the incompleteness of the Loveday et al. sample should allow to fully assess the consistency 
between the luminosity functions which can be derived from the two samples. 
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Figure Captions 

Figure 1): 

Confidence ellipses at la level for the parameters listed in Table 1, referred to the complete 
and the three incomplete samples of case 1). The parameters derived for the incomplete 
samples are clearly not consistent with the real ones. 
Figure 2): 

Parameter pairs (M* — a) found for each incomplete catalogue of case 2): open and solid 
circles refer to parameters derived with eq. (3) and (5) respectively. There is a clear separation 
between the two sets of points, being the solid circles very well consistent with the input 
parameters (star). The star (and the dashed ellipse) and the cross (and the dotted ellipse) 
refer to the parameters derived by Efstathiou et al. and Loveday et al., respectively. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/astro-ph/9403063vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/astro-ph/9403063vl 



This figure "figl-3.png" is available in "png" format from: 
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